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I. ABSTRACT 


Sh accurate method of calcwWlating the statistical monthly and anmaal * 
distribution of ambient temperatures on (treet Circle routes and over 

geographical regions is presented. An _ to operate computer program 
gives data in convenient graphical and tabular form for airplane system 


trade studies and performance calculations. ( ) ae 
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II. SUMMARY 


Accurate information concerning the statistical distribution of enroute 
ambient temperatures at airplane cruise altitudes is required for many 
trade studies and performance calculations, and is particularly important 
; when system performence is highly temperature dependent as, for example, 
. in airplene air-conditioning system studies. This document Zescribes 
a method of calculating the statistical monthly (January, April, July 
and October). and ennual temperature distribution on any Great Circle 
route for pressures of 300, 200, 150,°100 and 50 millibars, corresponding 
approximately to altitudes of 30,000 feet, 40,000 feet, 45.000 feet, 
53,000 feet, and @, 000 feet. 


The route temperature distributions are generated by a computer program. 
Inputs to the program are the coordinates of the terminal points. 
Output is in several forms, nanely: 

1. The route temperature distribution, showing the probability of 
exceeding any temperature, given in graphical form with the option 
of tabular form as well. 

The percentage of the total time that the temperature lies within 
discrete intervals of specified mediem and width, given in graphical 
form with the option of tabular form as well. 

The mean and standard deviation of the normal curve which neat. 
approximates the actual temperature distribution, and the error 
associated with the normal curve approximation. 

The graphical output is illustrated in Figures 6 and 7 for the 

Johannesburg to London route at 30.000 feet. 


Meteorological data based on records compiled over long periods for a 
selected global network of points furnish the basis for the program. 


Sections of an existing program (Boeing Document D6-6833TH, Program 
No. TAPOQ3) are employed, as a subroutine, for the determination of mean 
temperatures and standerd deviations at equidistant points 100-200 ' 
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nautical miles apart along a Great Circle route. Tr? results were ©. 
then used in the earlier program to generate the route mean terperature 
and the mean value of the enroute standard feviations for meny routes, 
as described in documents D6-15650-2 and D§-15650-4. . Presentatiun of 
results in this form has the disadvantage that the averaging process “ 
obscures the effect of temperature variations along a route. Furthermore, 
representation of the temperature distribution by @ normal curve results 2 
in a loss in accuracy. For some cases, and particulerly for short 
East-West routes and for all routes at 40,000 feet, enroute temperature 
variations are small and the route-mean results will closely approximate 
tha results developed in this document. However in many cases enroute - : 
temperature variations are significant, and misleading information about 
temperature extremes will be obtained if the route-mean resuits are 
employed to describe the enroute temperature distribution. 


Temperature distributions for specific routes may prove unnecessary 
voluminous to a user who only requires information for a specific region 
of the earth. Therefore results are presented for four geographical 
regions of interest. namely: ; 
1. The continental United States. 
2. The continental United States combined with the North Atlantic 

and Central and Southern Europe. 
3. The polar region of the Northern Hemisphere. 
4. The equatorial region. 


Temperature dietributions for these regions are compared with the | 
U. §. Standard Atmospheres, and it is shown that the Hot and Cold 
Atmespheres are encountered only over limited and often remotely 
located areas. 


The advice and cooperation of Mr. Neal M. Barr (Atmospheric Engineering 
Group. Aerodynamic Staff) is ecknovledged. 
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III. INTRODUCTION 


The embient air temperature influences the performance of many sirplene 
systems. Accurate information about this temperature is, therefore, of 
prime importance not only for trade studies but also for performance 
calculations. 


Many different methods of estimating representative ambient temperature 
distributions have been used. One approach is to consider just the 
defined values of Stantard, Hot and Cold Atmospheres, and to assume 

that these occur in fixed proportions (90 percent, 10 percent, 0 percent, 
for exemple). Such an approach ie clearly very arbitrary and is unlikely 
to yield accurate results. Some of the ahortcomings of this method are . 
avoided ‘by using the route-mean temperature data which have been compiled 
for weny Great Circle routes (Reference 1). Route data in Reference 1 
are given in the form of the route mean temperature, the mean value 

of the enroute standard deviations, and the temperatures that are not 
exceeded 50 percent, 75 percent. and 85 percent of the time. ‘The effect 
of inflight temperature variations is not reflected in these data and ‘ 
therefore information about temperature extremes calculated from the given 
mean temperature and standard deviation values, may be misleading. 


The method of calculation described in this dovument avoids the 
approximations introduced in employing the route-mean data of “¢ 
Reference 1 to determine the enroute temperature distribution of a 
studied route, and the error involved !n approximating the actual 
distribution by a normal curve. Both the developed method end the route- 
mean method use the same meteorological data. The deviation in the information 
about temperature extremes between the two methods is dependent on the 
magnitude of the enroute temperature variation, and this is inturna _ 
complex. function of the route, the altitude and the tine of the year. 

In general, the difference increases with length of route, is greater 
for Sorth-8outh than East-West routes, and is comparatively high at an 
altitude of 53,000 feet and low at 40,000 feat. 
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IV. CALCULATION OF TEMPERATURE DISTRITWTION ON A GREAT CIRCLE ROUTE } ee 


A. MBTEOROLOGICAL DATA 
Meteorological data, based on recordé complied over long periods, have been 
obtained for a network of 1117 points covering the surface of the earth. 

The points are located at every latitude which is en exact multiple of five ’ 
@egrees. Within 60° latitude of the equator, the points are locsted 

at every longitude which is an exact multiple of ten degrees; on those: 
latitudes which are farther than 60° from the equator the points are 
located at every longitude which is an exact multiple of twenty degrees. . 
Bach pole is represented by one point. ‘The mean and stander deviation of 

@ normal distribution fitted to the actual temperatures over a period of 

a month ere recorded for each point. Date are available for January, 

April, July and October, these months being assumed to be representative 

of the seasons; and for pressures of 300, 200, 150, 100 and 50 millibar . 
corresponding approximately to altitudes of 30,000 feet. 40,000 feet, - | 
45,000 feet. 53,000 feet end 68,000 feet. Thus a total of kk ,680 deta ig 
values are available. 


B. GREAT CIRCLE ROUZE CALCULATION a 
If the geographical ccordinates of the terminals of a route are 
specified, the Great Circle route may be determined by atenderd 
methods (Reference 2 for exazple). ‘The coordinates of most major 
airports can be found in Reference 1. 
C. CALCULATION OF ROUTE TEMPERATURE DISTRIBUTION 
' If a route is divided into 4 number of equidistant points, the mean 
temperature (m,), and standard deviation (0; ) at eech point (for e 
‘given month and altitude) may be obtained by interpeleticn from 
edjacent data points (Reference 2). Then if the tempersture ct any 


mM; ", 
m, 0; oo 
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point (4) is distributed normally, the probability that the tempers- 
ture T at this point is less than or equal to any specified ten- 
perature W is, by definition: 


. B ~m\™ 
PUT 1) =F, (WV), =a GS) at (1). 
a (VB eh 71 


a 
The probability that the temperature #% is exceeded is given by 


P(T >v)=1-F,, (OW) (2) 


If the distance between points is chosen so that the temperature 
and standard deviation can be assumed constant during each flight 
increment, then the probability for the whole route of any specified 
temperature 7% being exceeded is given by 


Pe (T> 29) = 1- Fe, (18) = qlee cee 5) 


An example of this method, using hand calculations to determine 
the cumulative probabilities, is show in figuree 1 and 2. In 
order to limit the number of calculations, the example uses larger 
temperature steps and fewer significant figures. The difference 
between the temperature distribution at selected points and the 


route temperature distribution jis itlustrated in figure 2. 


The probability of the temperature falling within a specified 
interval is given by: 


P(YeT < VtAV) © Fe (VtAv) - Fe (4) © 
Assuming that the months of January, April, July and October are 
representative of the four seasons, the annual probabilities may 
be obtained. Thus: ss 
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2 (1-629) 
Pan (T 2B)! It Fea (Be 


Pan(U<T <VtAV) = Fay (W+A)-Fe, (2) io 


D. FITTING OF A NORMAL CURVE 
For most routes the temperature distribution closely approximates 
a normal curve. Since normal curves may be defined by* specifying 
just tne-mean temperature and standard deviation, it is worthwhile 
determining the normal curve most closely approximating the actual 
distribution. 


The mean temperature of the normal curve is calculated from the 
ac'sl temperature which is exceeded 50% of the time and from 
temperatures having the defined probabilities of f c/z and Lo, 


thus averaging out discrepancies in any single value and giving: 


Meee ( Ue t Beg t Yen + Von, + U,) (7) 


The standard deviation is found by a similar method, averaging 
values at tO/2,4 O° 4.......1% 30, and giving: 


z— V-m + - 
Oe = 7g] DCO, /2) + (CH -mel/2) (8) 
? 
The mean and maximum deviations in (1—-Fa(W)) between the actual 
distribution and the fitted curve are calculated in order to indicate 


the accuracy of the curve-fitting. 


E. CALCUIATION OF A PROBABLE REGIONAL TEMPERATURE DISTRIBUTION 
The available meteorological data can also be utilized to generate 
a probable temperature distribution that would apply to a geographical 
region of the earth rather than individual routes. “he less 
voluminous results applicable to an airplane typically operating 


over the specified region are obtained at the expense of averaging © 
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_ generate the regional temperature distribution based on the 
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the temperature profiles over the region. 


A method of calculation similar to Paragraph IV, C is used to 


Meteorological data points located within the specified region. 
The: spheroid shape of the earth is accounted for by giving » 


area-balanced weights to individual points. Results are given in 
Appendix A. 
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DESCRIPTION OF PROGRAM 


A complete listing of the computer program for calculating enroute tem- 


perature distributions is given in Appendix B. it can be seen from the 


listing that the program is divided into a main program, five subroutines 


and & function subprogram. The main program handles input and output 


routines, *nd directs control to the subroutines and the function sub- 


program. The meteorological data is stored on a tape-unit (called TAPE 


for cross-reference) assigned to the author. The subprogrems have the 


following functions: 

1. Subroutine TEMP. 
Sections of an existing program (ref. 2) are used for generating 
mean temperatures and standard deviations at equidistant points 
100-200 nautical miles apart along the specified route. 

2. Subroutine DISTR. } 
Determines the statistical temperature distribution for the whole 
route from the mean and standard deviation date for the individual 
‘points. 

3 Subroutine FITNC. 

Determines: the mean and standard deviation of the normal curve most 
closely approximating the calculated distribution curve. 

4, Subroutine PINT. 

Determines the probability of the temperature being within intervals 
of given median and width. The system of median temperatures are 
centered around the Standard Atmosphere temperature. 

3« Subroutine WBIN 
Prepared calculated data for input to the CDC-6600 system : 
tape-writing subroutine WRTETP (ref. 4). 

6. Function subprogram CNORML e 
This subprogram calculates the cumulative probability for a normal 
distribution, given tre mean and standard deviation and the region 

for which the probability is required (ref. document PS-497, Program 

Library Department). 
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Input for the program is dencribed in Appendix C, and sample data for the 


route Johannesburg to London are included for illustration. It should be 


noted that exactly antipodel . points must not be specified as terminals 


of a route since such points do not give a uniquely defined Great Circle 
route, Furthermore, meteorological data are not aveilable on the SO millibar 
(68,000' altitude) level south of latitude 55 8, A diagnostic ic printed if ' 


the Great Circle route enters this region. 


have been written so that not more than two altitudes can be handled in 
any one run, and routes are restricted to a maximum length of 11,000 
nautical miles. A 
Output data are in the form of graphs and tables and are self-explanatory. 
The plots are not labeled by route, but may be identified by the fact 
that they .are generated in the same order that the routes are inputed. 
Tabular and graphical output corresponding to the input of Appendix C are 
shown in figures 5 through 7. 


Appendix D shows a complete deck assembly and lists the required control 
cards. The Fortran nomenclature is shown in Appendix E. The core storage 
requirement is 110000. and the central processor time is approximately 
equal to 


In order to reduce core storage requirements the DIMENSION statements 
CPT = 21.5 + 3.8 x NALT x NROUTE seconds 
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VI. EXAMPLE 


ROUTE TEMPERATURE DISTRIBUTION 

Use of the program may be illustrated by considering a specific route, 
and Johannesburg to London at 30,000' has been selected. As this is 
a long route (4896 N.M.) and in a South-North direction, considerable 
enroute temperature variations can be expected. This is confirmed by 
reference to figure %, which shows isotherms at 30,000’ in January. 


Figure 4 shows that the standard deviations varies noticeably too. 


The mean temperature and standard deviation for 26 equispaced points 
along the route are shown in figure 5, which also shows temperatures 
having » 50, 75 and 85% probability of not being exceeded. These 
values could also be obtained from the graphical plot of figure 6, 
snd can be compared with values obtained using the route-mean method 


of reference 1. For the month of January, the comparison shows! 
‘ - 


Prebability of penroute Route-Mean 
eaperature 
Not Exceeding Method 
50% -32.3F 
75% ~24 OF 
85% -22.6F 
' 


Figure 5 shows the mean and standard deviation of che normal distri- 
butio™ curve most closely resembling the actual temperature distribu- 
tion, and also the error associated with the normal curve approximation. 
It can be seen that in this case, the normal curve is not a particu- 
larly good approximation. A comparison for the month of January with 
the values Jisted in plete ebiich Sel Phd 


Temperature Route-Mean 
Method Method 


Mean Temperature -76.1F -38¢ 
Standard Deviation ~13.4F 
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It can be seen that the route-mean method gives misleading information 
about the extremes of temperature that will be encountered on this 

route and this method, originally not intended to describe the enroute 
temperature distribution but frequently used in this capacity, should not’ 
be used where considerable enroute temperature variations can be expected. 
It is interesting to apply the different methods of determining route 
temperaiures to the calculation of the operating penalty of a systen. 
Consider, as an exemple, the ram drag of the model 747 air conditioning 
system at .85 MN cruise. The ram drag is plotted in Figure 8 as a function 
of embient temperature, with temperature histograms based on the 
statistical and route-mean methods superposed. The average route ram drag 
may be compared, giving: , 


» ‘Average Drag 
Method ; - 1b/Airplane ; 
Enroute temperature 799 
Route-Mean ; 396 
Standard Day 90% 
Hot Day 10% M5 


It can be seen that in this case the route-mean and Standard/Hot day 
methods seriously underestimate the average drag, and hence the 
average airplane penalty. 


REGIOHAL TEMPERATURE DISTRIBUTION 

Appendix A introduces the use of the available weteorological glata to 
generate a probable temperature distribution applicable to a specified 
region of the earth, thereby eliminating the high degree of idealization 
inherent to the use of standardized atmospheres. To illustrate the use 
of the generated data, an example involving the identical airplane- 
system as in the previous example will be presented. Two of the selected 
geographical regions in Appendix A (Continental United States and 
Equatorial region) at the 30,000 foot altitude will be studied. 


From the annual temperature distribution curves for the two regions 
at the given altitude (Figures 14 and 17) one can easily obtain the 
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probabilities on 5 degrees Fahrenheit intervals for calewlétion ” 

of the average ram drag, as was tone for an individual youte - in 

the previous section. Figure 10 gives the probability histograme 
for the two regions superposed on the ran drag curve. Cmmnrlael ot 
the average ram drag gives: 


Average drag 
1b/airplane _ 


Standard Day 90% kis 
Hot Day 1% 
Probable temp. distribution over 


the Continental United States nal 


Probable temp. distribution over 
the Equatorial region — 1260 


It is seen that for the given exemple, the adopted combinations of the 
‘standardized atmospheres yields an acceptable value of the average 
drag if applied over the Cantinental United States. If. ancther 
airplane system, having a different temperature dependence, was being 
‘studied, this might not have been the case. Taken over the 

equatorial region it is found that 90% Standard Day, 10% Hot Dey 
seriously under-estimates the average drag, in spite of the fact 

that Hot Day (Hot Day = Std. Day + hOF at 30,000 ft.) 16 found 

"to be exceeded < 1% of the year when focusing the entire equatorial 
region of the earth. Likewise it can be seen that temperatures as 

low as Standard Day are extremely unlikely to be encountered over 

this area, making an assumption of 90% of the time at Standard Day 
unrealistic. Fanally if the interest is founded on locelized areas 
where the concept of Hot Dey applies (areas recording teseratures 
exceeding Hot Day 10% of the time), one will find that such as 

area at 30,000 feet exists over South East Asia. Calculations of the 
saverage drag based on the probable temperature distribution over this 
area will result in a value exceeding the one obtained for the Equatorial . 
region, showilg that if 10% Hot Day io a realistic figure, one can not 
designate Stamfard Dey condition to the remaining pontion of the tine. 
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VII. ComCTas ams AND RECOMMENDATIONS 


” The method for ‘determining enroute semparetnres developed: in this 


docunent genérates results which are limited in eccuracy only by. 
the accuracy of the meceorological data.on which they are based. 
The computer progrea is simple to use end presents results in 
convenient tabuler and graphical forms. Therefore, it is recommended 
that the statistical method should always be used for generating enroute 
temperature distributions. For some routes, and perticularly for 
short East-West routes, the route-mean method (Reference 1) gives 
acceptable approximations of the enroute temperature distribution. 

' 
Use of the regional temperature date to designate a probable temperature 
distribution eliminates the high degree of idealization inherent to the 
use of standardize’ atmospheres, thereby greatly imprewing the accuracy 
end reliability of employed data. It is recommended that the regional 
temperature data be used instead of standardized atmospheres in . 
determining typical operating temperatures for airplanes ani oe tac 
systems. 
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THE NORMAL DISTRIBUTION CURVE FOR A NUMBER OF 
SELECTED POINTS IS SHOWN. BELOW. THE CURVE 
CNOAI- NORMAL ) BASED ON THE CALCULATED DATA JA) 
FIGURE. GIVES THE DISTRIBUTION FOR THE EASTIZE 
ROUTE. pat 7 


_PEOBABILITY To EXCEED 


© ~50 -45 -40 -38 -50 -25 
TEMPERATURE ~ °c 
£ op) san 
EXAMPLE: CP.) J 2 a 


e 
“2-45 2b 


TEMPELATUGCE OISTEIBUTION 


aera | 
areo | CuRVeS 
ae 


TO O73 B-N5 


A 


~—— 


AVERAGE TEMPERATURE AT 300 MB. IN JANUARY 
UPPER AIR TEMPERATURE OVER THE WORLD 
jie N. GOLDIE, IG. MOORE AND &.B. AUSTIN 


+ 


| 


iy 
; 
; 
me Ie 
—— 
™ . 
. 
ai eee 
- . t 
mee’ 180 
tae we ary 
J te 
a 
; 
: 


oa 


t 


* 
4 = 
r 
ry ’ } 
’ < ‘ 
’ , 
; - 
n Ps 
' ' 
. s ,?. . 
. i . 4 iP 
r nal % —_ a) 
Ty EE AR sf | | 
“Z iy 4 
J - = ss 
os - 7 § - 
. 
A ‘ 


—~---+—- — — 
——e % 


“—-. 
“ 


—<- oo Or 
‘ 
t 


}? 


STANDARD DEVIATION OF TEMPERATURE AT 300 MB. IN JANUAI 
UPPER AIR TEMPGRATURE OVER THE WoeLDd 
Nl. GOLME, JG. MOORE 4ND EE. aysnn 


Se re 


i 


! 


20 
‘ oa 7 
_ 
{ 


te 
’ 
\ 
Ny 


U 
; 
4 “a 
} | 
GED land over 1Q.000 Fr. 
A " Ti 140 l 


60 


i 
| 
| 


xy 


| 


| 


Case WA, 1 oul oF NunTie ' 


JOAN £81108 10 LANDON enGA N,4, 
POSdqeesorcoocsoccoosoccosooooaneoosecolecees 


STan?y POINT Latifwe © oft, ft 
LONGITUDE © oFR.1% 
and PulN? Latitude © 61.2% 
Lonettvog © off 


OLTstuuce econ"? 
PY 


oc Toute 


330,6/902 
© 39.0736 
039.%72—8 
oP PLL the 
oS2ohs2e4 
odbetseok 
odecfsket 
050.07200 
o84%e67) 048 
o2%.67108 
229.9 sh0? 
o30.bs1,7 
od Oshe? 
o43,071.A 
32 ,0/1.80 
033,0716% 
039,271.69 
oJ0,872,3 
o38,472,3 
o50,677,8 
©90,9/5.3 
2b017509 
00} ,07%,! 
©02,6/0.2 
@39,370,0° 
006,0/%,0 


00/51 0EV, POW POINTS ALONG fhe wQUTE. eee CELCINS eee 
OFST wm, ganuaty APOE JULY 
0 032,777 0% os: 97908 ©37,A73_% 
ee: 3397706 290.9779 29SeF Vel 
302 CPL YY 2a 8 092007207 ©33,A79," 
daa CP Uy Ye a8 | 9940172 e% ©33.07/20% 
7e3 38/200 233037204 233077723 
ory oINnersie¥ 2300778? CDN 
gf 030% sieo? 2390027260 92 0hs2 00 
137) 031 .O7h 06 24,771.69 o32.371.9 
yee? eJteSsiee o2vedsio? Poh sioA 
saved oP .8sieb ~2%eesie? De Psd eA 
1980 037 .9shee o2%,%7i06 adioesto? 
aise 033,% rice o24,0s/1.4 od ed shah 
2380 096,% sic ©30.37366 ©3947 718 
gdee 03%, 8710? o3) 2sie? 290037b09 
avez 03777108 32,0740? @29,972.0 
2930 ©39,87700 099 .A71 04 @29 7208 
3333 062,372.01 o39,h72.) 299,572.09 
332° 263,07 208 o3f W779 ©30,3/7.9 
9d2o 040,970! 040.279.5 @33.877.8 
37a: o6D,8 776% o62,674.9 0359.9 7900 
3937 oOo TsVo2? 360,973.29 239.9730) 
YT¥) art SPL YA Oh 709,07 306 ©30.07/3.? 
een 068, 37900 ©6%.1 7000 O37 277344 
ebee o0%,h ocd 000,000? ©34,873.% 
oreo 000,97 390% o0f,t75e% ©3¥,873.6 
onto 009.97 90% oe! D730? eel so? 
TEMPEMATUNE UST IN FiuLOcING Tawrs IS Coe banMRENNELT eee 


TEMP, (MAT CORUESPUND 10 GIVEN PHUN, OT TH KF FECEFNFO 


PHOW, VALUE Januaiy aPeQtL JULY 
038 032.3 ou, o2!,t 
fd 020,99 @22,9? #2), 
205 22.6 20,9 e2). 


a 
. 


FET a eOnmaL CURVE 10 ImE MOUTF NISTH, 


JANUGR? aPelL JULY 
ap any Temes @imebi 532,00 21,95 
ST. Otve & 13ed0 1b.70 9.3an 


Hest nance OF PITTED Cumvt fo ROUTE DISTH, 


ncineew 
o29,0 


o?),? 
o21,% 


AC TOFU 


e161 
9.41 


OCVialsOn IN PROBACIL ITY vaLUEs BUM CALC, TEMPs RATURES 


Aull. “cugan Ntv.s 056 2038 oA? 
Manin Uv, 8 ohd: elle 2954 
Figure 5 


0019 
ea 


ANNUAL 


33.40 
10,06 


e0ev 
2046 


F cee 
, sey hy 
: ee 


Re ‘2 


¥ 
> 


i. ‘ago 
we ft Pee 
Me hee * ’ 
iat ee: ce 


NG ay 
Tees 

Late 
ree re 


z 


‘ 


* 


. 
. 
. 


rs 


sy 


PROBABILITY. TO ExCEED 
r 


. 


& ve wh ay : 
. any rae 7 ‘ 
a) . i bs . 
wy . ‘ .- e 
F : * 
« . ‘ . “ —_- *- * 
aaa Ze Na 4 ree f 
s sy : 


rf 


--* +e 


SEASONAL ANG ANNUAL OISTR, 
FOR GIVEN ALT. AND ROUTE 


TIL aU ae 


tee tae 08, ess Pre ee 


‘WOHANNESBUEL To LONDON 
S000 FT ALTITUDE 


Phat a 


"3 ~ 
Keio ee is 
aerate. 
Tage a0" 

. we 

sa: tle 

Hdl «a 

eed. hey ve 


Eee 
Kips new 
ane 
tof SAS 
Mad ‘ 
his . ry 
phe ie 


; 
‘PROBABILITY ON INTERVAL 


NS + ‘ 
cme 
‘et 
y Reuw 6 K 
a atts “ 
Tee 
ae on 
b as é 
% 
Agee ; 
tee " : 
34, 
Lay 
e ’ 
‘ 4 
i . 
° 
t 
‘ 
§ 
‘ 
si 


: TEMP INTERVAL REF. TO STO. 


PROBABILITY ON INTERVALS 
OF “SPECIFIED \ 0TH 


vue MPLIELIVAS . 


ia Bement qe are ee a ae 
peepiaeat HHH f fare 5 le i i eaii Hi ait por rare | ea 
A oe jon Pa ie si SERRE vet el eae ee PR ete 

ial aa S A se ain if: tT ae eeu gee 

a ggee oe Se 

fo sree 


ia Bey a rc a 
0 ee en ea 
ara ee eee eee 
ae 
ane eee 
L ete 
Soe: ea 
ae a ie Li. 


oe ae wae aoe 
as a, ee oo 
Hu ag 


ae _ 
He af Ge a. oo E 
a a cope oe cue 
é SSHFSEE iad itu fae SPR aT HGE Gh 


eG 
aye a ii . 


Hi He 
Ht 
eae 


lg 
ane 


THE BOEING COMPANY 


TO 41 C-R4 KeE Pyare 


ae 4 


aabeete 


ce a: nee gage 
ea 7 
Hi Hatt dip 
ee ee eee Oana 
Lt [eee Tae Sete Hane 
2 Ee a ae 
Sel one EH ee 
eee eneee lat 


fo ar A 

| Re Ei a fp eeiaes att ae 
He CU 
Bhi esl EB TESS He Pope a RAED rete a 
7 PEAR Tpate yy aa i a aire 
ie pprpd bs ee ‘ 
‘RS AR Ee ee aT 
7 ete elo g eh tit eA 


bee 
rthigsy 


ti ses if a 


t. i PG ME 

a Se age a HL aa A A 

4 a ee 

f H iF ER Fe EE fe ms ane “FREER TER SR a es Heed pitt fli 

fa ce EERE 
ee ee 


aia i ag Ps 
ae a ae : 


ae 
oe 
i He eee 


at Hie 


% S Be i BREE 

SiH SESH ERE HEE Sat ad a si — fH atte 
ae ’ ce 7 a 
Ha ae : oe 

ASR ES ooo PER oe ses eb SeOAOaE HR Se VP ef ae 
T41 AI CONDITIONIAIG PACK RAM Ale 
ORAG, WITH REGIONAL TEMPERATURE 
DISTRIBUTION) D6 -58402 | 


THE BOEING COMPANY 


TO 4 KeE eee ‘eee. 


IX. REFERENCES 


1. Document D6-15650-2, World Enroute Temperatures, N. M. Barr,, 
E. M. Hansen, March 1967. 


2, Document D6-6833TN, Enroute Temperatures, L. A. Rasmussen, Dec. 1963. 


3, Mathematical Handbook for Scientists and Engineers, Korn and Korn, 


Mc-Graw Hill Company, AIR 005471. 
4, Document D6-19175TN, SC-4020 Plot Package, David D. Renhed, Dec. 1965. 
5. U. S. Standard Atmosphere, 1962. 
6. MIL. STD -P1lOA. Climatic Extremes for Military Equipment. 


7» Atmospheric Environmentai Test and Design Criteria, Diamond M. and 
Essenwenger, O. M., Tech. Report SELWC-M-1, Sept. 1962, Army Electronics 


R&D Activity, White Sands Missile Range, N. M. 


AD 1846 0 


REV SYM . BOEING \~0. Do-58402 > 
PAGE 29 
6-7000 


AO 1846 © 


APPENDIX A 
USE OF THE AVAILABLE METEOROLOGICAL DATA 
TO GENBRATE REGIONAL TENPERATURE DISTICIBUTIONS 


The concepts of Standard, Hot and Cold Atmosphere (References 5 through 7) 
are commonly coabined to designate a probable temperature distribution 

for trade studies or performance calculations on the typical operation 
condition of an airplane system being studied. Any such combination of the - 
standardised atmospheres (1.e., 90 percent of the time at Standard Dey and 
remaining 10 percent at Hot Day) is likely to involve a substantial degree 

of idealization of any probable temperature distribution and the risk for 
tisinterpretations becomes apparent. If Standard Day condition is encountered 
during the major part of the time, ar in above exemple, the probability 

of experiencing Hot Day is extremely remote. On the other hand, if’ the 
interest is focused on the localized areas where Hot Dey temperature actually 
is being exceeded the specified percentage of the time, it is improbable 

thet temperatures as low as Standard Day will be encountered. 


Using the available meteorological data, it is possible to generate a 
more realistic temperature distribution curve applicable to a specified 
region of the earth, thereby eliminating the uncertainties inherent in the 
use of the standardized atmospheres. ‘The developed data becomes applicable 
to an eirplane in traffic of about equal density over the specified region. 
Date for four regions of common interest has been generated. These are: 
1. Continental United States with adjacent border areas (52.5N to 22.55 
latitude and 65W to 125 W longitude). 
2. Area 1 combined with North Atlantic, Central and Southern Burope 
(Area 1 plus 62.5 to 37.5N latitude and 258 to 65W longituie). 
3. Polar region of Northern Hemisphere (90M to 62.5N latitude). 
4. Equatorial region (22.5N to 22.58 latitude). 


Data pertinent to any other region of the earth of interest can, upon 
request, be generated. % 
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Pigures 10 through 13 give the temperatures that are. seeks to ve ees 
exceeded 99, 90, 59, 10 and ) percent of the year over the different 9. 
specified regions. The three standardized atmospheres (Standard, Kot and: gota” 
are superposed for comparison and reference. The temperatures at the =. is 
different levels where calculetions have been carried out (10,000, 28,000,” 
30,000, 40,000, 45,000, 53,000 and 68,000 feet) are interconnected in 

the most probable way. 


Figures 14 through 17 give the complete seasonal and annue) (arithmetic 
mean of the four seasons) temperature distribution at 30,000 ft. altitude. 
Similer curves can be made available for any region and altitude. ‘ 
Figure 18 is a summary table of data from Figures 10 through 13, 
everything being referred to Standard Day condition at each altitude. 
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“v APPENDIX 


INPUT DATA DESCRIPTION AND ARRANGEMENT 


The table belew describes the preparation of the data inyzut carde. 
LOCATION INFUT MODE OF 


CARD COLUMN ITEM” NUMBER DESCRIPTION : " 
1 1-6 ‘NALT INT. Number of altitudes tb be 
; used throughout the caleu- 
lations. 
7-12 ALTI(1) INT. First altitude (feet). for 
which data is to be caleu- 
lated. r 
13-18 ALTI(2) | INT. Seconé altitude in caleu- 
lation. ALTI(2) mut be. 
: larger than ALTI(1) (altt-"° 


tudes entered in order of | 
increasing numbers). AL®T(2) ) 
omitted if NALT equal to one. 


2 1-6 °  OPT(1) —_sINT. Option; = 1 if printing of 
all caleulated data, = O if no 
7-12 OPT(2) INT. Option; = 1 if temp. output in 
, Fahrenheit, « 2 if Celeiua 
13-18 OPT( 3) INT. Option; = 1 if probability on 
intervals desired, » 0 if not 
19-30 TINT REAL Width of temp. intervals, can 


be omitted if OPI(3) = 0, 
OPT({2) controls input unit. 


3 1-6 NROUTE INT. Number of rovtes to be calcu- 
lated in the vpecified sanner. 


4 1-80 COMMENT ALPHAMUM. Description of route; This 
card must neat be omitted. 


5 1-12 GIAT(1) REAL Start point latitude in degree 
and minutes to the conven~ 
A tional notation; north lati-. 
tudes positive,south negative. 


13-24 GLIN@(1) = REAL Start points longitude in de- 
grees and numbers to the con- 
ventional notation; west long- 
itudes positive, east negatives 


25-36 GIAT(2) REAL End point latitude 
- 37-48 GING(2) = REAL End point longitude 
3 
oO 
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Input data to process first work is now complete. For additional 


routes, repeat cards + and 5 until specified NROUTE are imputed. 


Input data to generate the sample case output shown in figures 5 


through 7 are listed below: 
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Johannesburg to London 4896 N.M. 
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‘The following table gives the nomenclature for the main program and 


SYMBOL DESCRIPTION UNITS 
A (1I,J,K) Meteorological data. Coded in the C (Celius) 


APPENDIX. E 


FORTRAN NOMENCLATURE 


the subprograms. Only the more significant terse’ are included. ‘he 
nomenclature for terns, whose values are stored ih the COMMON région 
of the computer, are given first. In the table, dummy subscripte z, 


J and K are used to denote 1, 2 and 3rd dimension respectively. 


NOMENCLATURE FOR THE COMMON STORAGE 
FORTRAN ¢ 


form taaabb, where taza is mean-temp. 
in tenths. and pb is standard devia- 
tion in .tenthe also. 

I denotes location on the earth 

J denotes month 

K denotes altitude 


DHBR(I) Temperatures calculated’ and used for C 
one altitude 


DIST(1) Cumulative distance from start’terminal N.M. 
to point i 


EAVR(I) Average error in fitted normal curve --- 


EMAX(I) Maximum error in fitted normal curve --- 
e 


ITHI Upper boundary of temp. range con- C 
sidered for one altitude. 


ITO Lower boundary of temp. range con- c 
sidered for one altitude. 


ITR Nuxher of temperature steps used in 
calculation 
I55S Indicator for region of earth where 


meteorologicul raw data are nct 
available (68000 ft, south of 55 S.lat.) 


Altitude 
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AA 


DESCRIPTION: 


Number of points i on the route 


Months number. 1 through 5 per 
following denotation. 


1 for January 
, 2 for April 

3 for July 

4 for October 
5 for Annual 


NINT 


Number of intervals employed for 
probability on intervals 


PINV(I,J) Probability on individual temp. N (Normalized) 
intervals of given width and 
median temperature. i 


PRI(I,J) Cumulative probability value at lower N 
boundary of interval 


' 


-PROB(I,J) Cumulative probability to exceed N 
temp. values given by DEGR(1) 


SCAILC(I) Standard deviation for fitted normal C 
distribution curve 


STDN(I,J,K) Standard deviation for the points Cc 
along the route. 


TIM(I) | Interval median temperature C 


TMN(1) Mean values for fitted normal C 
distribution curves 


TNMP(T ,J,K) Mean temp. for the points along Cc 
the route 

é 
TPR(I,J) Temperatures corresponding to .50, c 


-75 and .85 probability not to exceed 


USSA(1) United States Standard Atmosphere, 1962 Cc 


ADDITIONAL NOMENCLATURE FOR THE MAIN PROGRAM AND SUBPROGRAMS 
For complete nomenclature to subroutine TEMP, see Document D6-6833TN. 


Altitudes on which calculations are to Feet 
be carried out. 


ALTI (I) 
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FORTRAN 


CCF (Dummy ) 


OPpr(I) 
PIOD(I) 


PROD(I) 


| 

| 

| 

: 7™1(I,J) 
YSDL(I,J) 
| 


REV SYM 


Arithmetic stetenent functions for =; 2 

temperature conversion from Celcius 

to Fahrenheit. 

Latitudes of terminal points °.3 
V : 

Longitudes of terminal points de 

Error indicator for writing on o-- 

binary tapes 

Level (altitude) FEET 

Mumber of altitudes used in th “<= 

calculations : : 

Mmber of routes one 

Binary tape number --- 

Options for calculations oo- 

Arrey PIRV(I,J) tranformed to l= 5 

dimension 

Array PROB(I,J) transformed to °* | 

1=dimension 

Array THMP(I,J,K) transformed for 2=dim ¢ 

Array STDN(I,J,K) transformed to Cc 

2edin. 


BOEING | 10. 16-58402 
PAGE T1 


> 


6-7000 


. LI£? OF NOMENCLATURE z 
- Bymbol Description 
Fru (WH) Cumulative propability that T shell be less than ! 
or equal to 12 , for point "1° of a route. 
Fe (28) ' ' Cumulative probability ( T< WY ) for the entire 
route. 
Fra (Y): Fa(¥) for the annual case 
Ms Mean value for fitted normal curve. 
mi Meen value for point "i" of a route. t | 
n Number of points on a route, terminals included. | 
on] Standard deviation for fitter normal curve. | 
or: Standard deviation for point "i" of a route. 
Oso Temperature for Fp equal to 50 percent. 
ty Interpolated temperature used when fitting a normal | 
curve. | 
T Function that describes the occurrence of different ! 


temperatures UW at a given point. In this case 
a Gaussian function. 


vw Independent variable (temperature) 


AD 1846 0 


REV SYM BOEING \\1. %-58h02 ae 
PAGE T2 
6-7000 


AD 1846 0 


REV SYM 


Number 


10 
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12 


14 
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LIST OF ILLUSTRATIONS 


Title 
Illustration of Calculation Method, Part 1. 
Tllustration of Calculation Method, Part 2, 
ieee Temperature at 300 mB in January. 


Standard Deviation of Temperature at 700mB in 
Janvorye 


Program Printout (Sample case) 


Seasonal and Annual Distribution for Given 
Altitude and Route (Sample case) 


Probability on Intervals of Specified Width 
(Sample Case) 


747 A/C Pack Ran Air Drag. Route Temperature 
Distribution. 


747 A/C Pack Ram Air Drag. Regional Tempera- 
ture Distribution. 


99, 90, 50, 10 and 1% Probability of Exceeding, 


’ 


Temperature. Continental United Statee. 


99, 90, 50, 10 andl Probability of Exceeding, 
Temperature. Continental United States com- 
bined with North Atlantic, Central and Southern 
Europe. 


99, 90, SK, 10 anc.% Probability of Exceeding, 
Temperature. Polar Region of Northern Hemis- 
phere. 


99, 90, 50, 10 andl Probability of Exceeding, 
Temperature. Equatorial Region. 


Seasonal and Annual Distribution for Specified 
Geographical Region.Continental United States 
at 30,000 ft. 


Seasonal and Annual Distribution for Specified 
Geographical Region. Continental United States 
combined with North Atlantic, Central and 
Southern Europe at 30,000 ft. 


2h 


25" 


27 


28 


32 


33 
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Title 


Seasonal and Annual Dietribution for Specified 
Geographical Region. Polar region of Northern 
Hemisphere at 30,000 ft. 


Seasonal and Annual Distribution for Specified 
Geographical Region. Equatorial Region at 
30,000 ft. 


18 Deviations from Standard Atmosphere. ; Lo 


Deck Set-Up 
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temperatures and standard deviations at equidistant points 100-200 


II. SUMMARY 


Accurate information concerning the statistical distribution of enroute 
ambient temperatures at airplane cruise altitudes is required for many 
trade studies and performance calculations, and is particularly important 
when system performance is highly temperature dependent as, for example, 
in airplane air-conditioning system studies. This document describes 

a method of calculating the statistical monthly (January, April, July 
and October) and annual temperature distribution on any Great Circle 
route for pressures of 700, 500, 300, 200, 15Q 100 and 50 millibars, 
corresponding approximately to altitudes of 10,000 feet, 18,000 feet, 
30,000 feet, 40,000 feet, 45,000 feet, 53,000 feet and 68,000 feet. 


The route temperature distributions are generated by a computer program. 
Inputs to the program are the coordinates of the terminal points. 
Output is in several forms, namely: 

l. The route temperature distribution, showing the probability of 
exceeding any temperature, given in graphical form with the option 
of tabular form as well. 

2. The percentage of the total time that the temperature lies within 
diecrete intervals of specified median and width, given in graphical 
form with the option of tabular form as well. 

3. The mean and standard deviation of the normal curve which best 
approximates the actual temperature distribution, and the error 


associated with the normal curve approximation. 
The graphical output is illustrated in Figures 6 and 7 for the 
Johannesburg to London route at 30,000 feet. 


Meteorological data based on records compiled over long periods for a 
selected global network of points furnish the basis for the program. 


Sections of an existing program (Boeing Document D6-6833TN, Program 
No. TAPOO3) are employed, as a subroutine, for the determination of mean 


Reproduced by the 
CLEARINGHOUSE 
for Federal Scientific & Technical 
Information Springfield Va. 22151 
Reva wamina |v. 6-58402 
. Pace ) 


cn 


AO 1546 D 


REV SYM BORING \wo. -58402 a 


A 


IV. CALCULATION OF TEMPERATURE DISTRIBUTION ON A GREAT CIRCLE ROUTE 


METEOROLOGICAL DATA 

Meteorological data, based on recordé complied over long periods, have been 
obtained for a network of 1117 points covering the surface of the earth. 
The points are located at every latitude which is an exact multiple of five 
degrees. Within 60° latitude of the equator, the points are located 

at every longitude which is an exact multiple of ten degrees; on those 
latitudes which are farther than 60° from the equator the points are 
located at every longitude which is an exact multiple of twenty degrees. 
Each pole is represented by one point. The mean and stendard deviation of 
a normal distribution fitted to the actual temperatures over a period of 

a month are recorded for each point. Data are ava‘lable for January, 

April, July and October, these months being aszumed to be representative 

of the seasons; and for pressures of 700, 500, 300, 200, 150, 100 and 50 
millibar corresponding approximately to altitudes of 10,000 feet, 18,000 feet, 
30,000 feet, 40,000 feet, 45,000 feet, 53,000 feet and 68,000 feet. Thus a 
total of 62,550 data values are available. 


GREAT CIRCLE ROUTE CALCULATION 

If the geographical coordinates of the terminals of a route are 
specified, the Great Circle route may be determined by standard 
methods (Reference 2 for example). The coordinates of most major 
airports can be found in Reference 1. 


CALCULATION OF ROUTE TEMPERATURE DISTRIBUTION 

If a route is divided into a number of equidistant points, the mean 
temperature (m,), and standard deviation (0; ) at each point (for a 
given month and altitude) may be obtained by interpoletion from 
adjacent data points (Reference 2). Then if the temperature at any 


mM; mm, 
att B 
co; i) 
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Input for the program is described in Appendix C, and sample data for the 
route Johannesburg to London are included for illustration. It should be 
noted that exactly antipodal points must not be specified as terminals 

of a route since such points do not give a uniquely defined Great Circle 
route. Furthermore, meteorological data are not available on the 700, 

500 or the 50 millibar (10,000', 18,000' or 38,000’ altitude) levels south of 
latitude 55 S. A diagnostic is printed if the Great Circle route enters 

this region on those altitudes. 


In order to reduce core storage requirements the DIMENSION statements 

have been written so that not more than two altitudes can be handled in 

any one run, and routes are restricted to a maximum length of 11,000 

nautical miles. | 


Output data are in the form of graphs and tables and are self-explanatory. 
{ The plots are not labeled by route, but may be identified by the fact 

| that they are generated in the same order that the routes are inputed. 

| Tabular and graphical output corresponding to the input of Appendix C are 
| shown in Figures 5 through 7. 


Appendix D shows a complete deck assembly and lists the required control 


cards. The Fortran nomenclature is shown in Appendix E. The core storage 
requirement is 110000, and the central processor time is approximately 
equal to 


CPT = 21.5 + 3.8 x NALT x NROUTE seconds 
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